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1 .  INTRODUCTION 


The  efforts  described  herein  are  part  of  a  project 
to  create,  maintain  and  utilize  an  extensive  neutral 
atmospheric  density  data  base  qenerated  from  four  AFGL 
satellite  accelerometer  experiments  -  three  on  the  NASA 
Atmosphere  Explorer  (AE)  series,  AE-C,  AE-D,  AE-E,  and 
one  on  the  Air  Force  S3-1  satellite.  The  density  data 
were  taken  over  a  three  year  period  from  December  1973 
to  November  1976.  The  accelerometer  data  processing  sys¬ 
tem,  the  analytical  techniques  used  to  extract  atmospheric 
density  from  the  total  accelerometer  output  signal,  and 
the  creation  of  a  density  data  base  have  been  described 
previously  (1,2,3). 

This  report  describes  the  finalization  of  an  extensive 
neutral  atmospheric  density  data  base  taken  from  the  four 
satellite  accelerometer  systems,  the  development  of  an  em¬ 
pirical  model  of  these  data,  and  other  aspects  of  data 
analysis  support  for  AFGL. 

Studies  by  AFGL  scientists  utilizing  these  results  are 
described  elsewhere  (e.g..  Ref.  (4)). 


1 . 1  AE/S3  Data  Base 

The  AE/S3  density  data  base  includes  data  from  the  ellip¬ 
tical  orbit  portions  of  the  respective  missions  of  the  AE-C, 
-D,  -E  and  S3-1  satellites.  Table  1  provides  orbital  char¬ 
acteristics  and  dates  of  elliptical  orbit  data  acquisition 
for  each  satellite.  The  altitudes  at  which  data  were  obtained 
are  shown  in  Figure  1. 
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1. 2  New  AE-C  Orbits 

More  than  one  thousand  (1000)  orbits  of  AE-C 
elliptical  orbit  data  not  previously  processed  became 
available  and  were  processed  under  the  present  contract. 
These  data  were  processed  using  the  MESA  data  reduction 
system  (previously  described  in  Reference  (1)),  and 
calculated  density  values  were  stored  in  MESA  GU  and 
UA  files. 

The  new  AE-C  density  data  base  was  then  developed 
from  the  new  UA  density  data  using  the  techniques  des¬ 
cribed  in  Reference  (3).  This  data  base  was  generated 
on  the  Sigma-9  computer  and  was  written  to  magnetic  tapes. 
These  tapes  were  shipped  to  AFGL  for  insertion  into  the 
AE/S3  data  base  on  the  AFGL  CDC  6600  computers. 

To  accomplish  this,  programs  were  developed  to  insert 
(i.e.,  merge)  the  new  data  into  the  existing  AE-C  data 
base  for  all  altitudes,  to  quality  check  the  resultant 
merged  data  base  for  duplicate  data,  and  finally  to  re¬ 
merge  this  data  base  with  the  remainder  of  the  AE/S3  data 
at  all  altitudes.  This  new  composite  data  base  makes  up 
the  present  AE/S3  density  data  base. 

1. 3  Data  Base  Properties 

Figures  2-6  are  histograms  of  the  data  distribution 
as  a  function  of  geomagnetic  latitude,  local  time,  Kp 
index,  solar  flux  and  satellite  altitude.  Figure  2  shows 
the  number  of  data  points  obtained  in  each  20  degree  geo¬ 
magnetic  latitude  band.  Figure  3  shows  the  distribution 
of  data  points  obtained  for  each  hour  of  local  time.  The 
distribution  of  data  with  respect  to  geomagnetic  activity 
is  given  in  Figure  4.  This  gives  the  number  of  measure¬ 
ments  during  which  the  three-hourly  Kp  index  (with  6  hr. 
lag)  fell  within  a  specific  one  unit  interval.  Tn  Figure  5 
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Figure  2.  Data  Distribution  as  a 
Function  of  Geomagnetic  Latitude 
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Figure  3.  Data  Distribution  as  a 
Function  of  Local  Time 
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Figure  4.  Data  Distribution  as  a 
Function  of  Geomagnetic  Activity 
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Figure  5.  Data  Distribution  as  a 
Function  of  Solar  Flux 


Figure  6.  Data  Distribution  as  a 
Function  of  Altitude 
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the  number  of  measurements  durinq  which  the  daily  solar 
flux  value  (with  one  day  laq)  was  within  specific  5  solar 
flux  unit  intervals  is  given.  Figure  6  illustrates  the 
number  of  points  in  each  5  km  altitude  file. 

An  example  of  the  data  base  is  qiven  by  showinq  results 
obtained  at  180  km.  Figure  7  shows  the  density  data, 
corresoondina  solar  flux  values,  ratio  of  measured  density 
to  the  J71  model  and  the  Kp  index  presented  as  a  function 
of  calendar  year  and  GMT. 


2 .  Model  Development 

2 . 1  J71,  J77  and  MSIS  Comparisons 

Histograms  of  the  measured  density  values  were  pro¬ 
vided  to  AFGL  to  permit  analysis  in  relation  to  three 
atmospheric  models  -  J71,  J77  and  MSIS.  Graphical  repre¬ 
sentations  of  frequency  distributions  of  percent  differ¬ 
ences  from  the  models  were  developed.  Figures  8-10 
display  these  histoqrams  for  the  entire  data  base  for  each 
atmospheric  model.  Figures  11  -  13  display  these  histo¬ 
grams  for  200  km  data  only.  Frequency  distributions  of 
the  data  were  described  in  terms  of  the  mean  value  and  the 
second,  third  and  fourth  moments  about  the  mean. 

2 . 2  Empirical  Model 

A  multiple  linear  regression  analysis  program  was 
utilized  (3)  to  support  the  AFGL  development  of  an  empir¬ 
ical  model  of  atmospheric  density  variations  in  terms  of 
latitude,  longitude,  local  time,  geomagnetic  activity, 
solar  flux,  season  and  semi-annual  effects.  This  program 


Figure  7.  Data  Base  at  180  km:  Solar  Flux,  Density 
Ratio  to  J71  Model  and  Kp  as  a  Function  of  Time 
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Figure  10.  Frequency  Distribution  of 
Percent  Differences  of  all  Data  from 
J77  Model 
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Figure  11.  Frequency  Distribution  of  Figure  12.  Frequency  Distribution  of 

Percent  Differences  of  200  km  Data  Percent  Differences  of  200  km  Data 

from  J71  Model  from  MSIS  Model 
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Figure  13.  Frequency  Distribution  of 
Percent  Differences  of  200  km  Data 
from  J77  Model 
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was  expanded  to  have  the  capability  to  include  a  much 
larger  data  base  -  specifically,  the  composite  AE-C,  AE-D, 
AE-E  and  S3-1  data  base  -  and  more  independent  variable 
terms.  We  then  further  expanded  this  program  to  have 
the  capability  to  perform  stepwise  linear  regression  anal¬ 
yses  on  the  data  base  to  determine  the  amount  of  variation 
in  the  data  which  is  described  by  each  term,  the  signifi¬ 
cance  of  each  term  and  the  order  which  each  term  enters 
the  model.  Under  the  present  conf iguration  the  multiple/ 
stepwise  regression  analysis  program:  (a)  has  the  capa¬ 
bility  to  process  a  data  base  at  least  ten  times  larger 
than  the  original  AE-C  data  base?  (b)  contains  over  one 
hundred  (100)  independent  variables  from  which  may  be 
chosen  the  variables  to  be  used  in  the  particular  model 
to  be  run?  and  (c)  allows  up  to  fifty  independent  vari¬ 
ables  per  model. 

An  expanded  version  of  the  empirical  model  being  de¬ 
veloped  in  support  of  AFGL  contains  21  independent  vari¬ 
ables  and  was  run  at  each  5  km  altitude  from  140  to  240  km. 
The  multiple  regression  program  was  first  run  to  calculate 
the  phases  of  the  trigonometric  terms,  and  then  the  step¬ 
wise  regression  was  run  to  determine  the  order  and  signi¬ 
ficance  of  each  term  in  the  overall  model.  Table  2  dis¬ 
plays  statistical  data  related  to  the  regression  analysis 
including  regression  coefficients,  stepwise  correlation 
coefficients,  and  the  significance  and  order  that  each 
term  entered  the  analysis.  Although  calculated  for  each 
5  km  interval,  these  values  are  qiven  in  20  km  intervals 
for  this  report. 
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Table  2.  Regression  Analysis  Results 
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2  *  3  Results 

Utilizing  the  regression  analysis  we  were  able  to 
determine  and  display  amplitude  and  chases  of  annual  and 
diurnal  variations  (and  their  harmonics).  Figures  14  - 
15  display  some  of  these  results  as  a  function  of  satel¬ 
lite  altitude. 

Variations  due  to  semi-annual  effects,  Ap  values,  and 
solar  flux  information  are  given  for  200  km  data  in 
Figures  16  -  18.  Each  figure  displays  data  for  latitudes 
of  0  and  60  degrees. 

Because  the  empirical  model  is  based  on  data  at  selec¬ 
ted  altitudes,  it  is  necessary  to  perform  some  form  of 
interpolation  or  fitting  to  obtain  a  continuous  relation 
between  density  and  altitude.  We  have  chosen  to  fix  (i.e., 
hold  constant)  all  the  parameters  in  the  model  except 
altitude.  This  generates  a  density  value  at  each  altitude, 
and  then  these  values  are  fit  with  a  function  of  the  form: 

d (h)  =  ( c o  +  Cih  +  c2h  +  c3h  )  e  ^a°  +  a]h^ 

where  d  is  the  density  at  an  altitude  h. 

Figure  19  displays  the  resultant  model  fit  as  a  function 
of  altitude  for  March  21  (day  80)  and  June  21  (day  172) . 

The  difference  between  the  two  curves  is  a  consequence  of 
the  change  in  day  of  year.  The  scale  height  H  (=  'Vaj) 
changes  from  25.45  km  for  day  80  to  24.95  km  for  day  172 
which  is  consistent  with  the  larger  ranqe  of  density  shown. 
Table  3  displays  the  model  coefficients  for  the  curves 
displayed  in  Figure  19. 
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Figure  15.  Amplitude  and  Phases  of  Diurnal  and  Terdiurnal  Variations 
of  the  Data  Base 


ALMTUDt  IPNl  ALT 


DENSITY  CX10-13)  GM/CC  OENSITY  (X10-13)  GM/CC 


L5 


200  KM  COMPOSITE  DOTH  LAT=0 


DAY 


200  KM  COMPOSITE  DATA  LAT=60 


DAY 


Figure  16.  Semi-Annual  Density  Variation 
for  each  model  at  200  km  lor  0  decrees 
and  60  decrees  latitude 
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Figure  17.  Geomagnetic  Activity  Variations 
of  Density  for  each  model  at  200  km  for 
0  degrees  and  60  degrees  latitude 
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Figure  18.  Solar  Flux  Variations 
of  Density  for  each  model  at 
200  km  for  0  degrees  and  60 
degrees  latitude 
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In  order  to  support  studies  of  density  variation  as  a 
function  of  both  latitude  and  day  of  year,  a  set  of  con¬ 
tour  plots  was  generated.  To  accomplish  this  all  parameters 
except  latitude  and  day  of  year  were  fixed.  Figures  20 
and  21  display  density  contours  as  a  function  of  latitude 
and  day  of  year  for  Ap  =  0  and  Ap  =  80,  respectively.  The 
density  values  are  shown  in  units  of  the  ratio  of  the  em¬ 
pirical  model  density  to  the  actual  average  annual  density 
for  that  altitude,  with  a  contour  interval  of  0.025  in  those 
units . 


3.  Other  Contract  Functions 


3 . 1  EUVS  Experiment  Support 

Utilizing  our  knowledge  (gained  in  supporting  the 
MESA  accelerometer  experiment)  of  the  AE  data  system,  the 
AE  Sigma-9  computer  system,  the  AE  Data  Management  Facility 
(DMF )  software,  and  the  general  concents  and  philosophy 
of  the  AE  Aeronomy  Team  of  scientists  we  were  required  to 
support  the  efforts  of  the  solar  ultraviolet  experiment 
(EUVS)  scientists  at  AFGL.  We  utilized  the  remote  terminal 
dedicated  to  EUVS  to  provide  interactive  on-line  support. 
This  included  maintenance  and  execution  of  interactive  and/ 
or  batch  programs;  copying  various  output  files  to  the 
remote  printers;  promotion  and  demotion  of  telemetry, 
orbit/attitude,  geophysical  unit  and  unified  abstract  files; 
and  general  account  and  file  maintenance  to  insure  smooth 
execution  of  standard  EUVS  processing.  In  addition,  manual 
plotting,  hand  calculations  and  file  maintenance  programs 
were  written  to  support  this  requirement. 

EUVS  geophysical  unit  data  wore  written  to  magnetic 
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tape  on  the  Sigma-9  computer,  and  these  tapes  were  then 
shipped  to  AFGL  to  be  read  by  the  CDC  6600  computer.  We 
evaluated  various  tape  formats  and  wrote  programs  to  deter¬ 
mine  the  readability  of  these  tapes  for  the  EUVS  scientists. 

3. 2  AE-C  Re-Entry  Data 

During  the  re-entry  phase  of  the  AE-C  mission  important 
accelerometer  density  data  were  gathered.  The  data  were 
taken  over  the  final  twelve  days  at  a  rate  of  one  orbit 
per  day,  alternating  spinning  and  despun  mode  data.  Seven 
spinning  and  six  despun  orbits  were  loqged. 

Of  the  seven  spinning  orbits  two  orbits  contained  less 
than  three  minutes  of  usable  data;  on  two  orbits  the  in¬ 
strument  was  turned  off;  we  were  able  to  process  data  from 
another  orbit;  and  two  orbits  were  successfully  processed 
to  density  values. 

Of  the  six  despun  orbits  one  orbit  contained  less  than 
three  minutes  of  usable  data,  and  the  other  five  orbits 
were  successfully  processed  to  density  values.  Instrument 
bias  values  were  calculated  from  previous  high  altitude 
circular  despun  data. 

These  re-entry  data  will  be  incorporated  into  the  AE/S3 
density  data  base. 

3 . 3  Circular  Orbit  Processing 

The  obiective  of  the  AE  mission  was  to  study  pheno¬ 
mena  in  the  atmosphere  at  altitudes  above  120  km.  This  was 
to  be  accomplished  in  two  orbital  phases:  elliptic  orbit 
and  circular  orbit.  After  the  initial  elliotic  orbit  phase 
was  completed,  AE-C  and  AE-E  were  circularized  for  the 
second  phase  of  their  mission. 

For  circular  orbit  spinning  data  our  processing  tech¬ 
niques  to  calculate  density  (Reference  (2))  were  for  the 


most  part  unchanged.  For  despun  circular  orbits,  since 
instrument  bias  values  are  required  and  the  satellite  alti¬ 
tude  was  below  the  "non-drag  regions"  used  in  the  elliptic 
orbit  phase,  an  alternate  method  was  required  to  determine 
bias  values.  We  utilized  the  techniques  described  in 
Reference  (5)  to  calculate  bias  temperature  coefficients 
for  each  satellite.  Circular  orbit  temperature  data  was 
then  used  to  calculate  instrument  bias  values,  and  filtering 
techniques  were  used  to  determine  density  values.  Circular 
orbit  density  data  are  stored  in  the  AE  unified  abstract 
files . 

4.  CONCLUSIONS 

An  extensive  data  base  has  been  developed  utilizing 
accelerometer  measurements  obtained  with  four  low  altitude 
satellites.  Utilizing  this  data  base  an  empirical  model 
of  neutral  atmospheric  density  has  been  developed  by  mul¬ 
tiple/stepwise  linear  regression  techniques.  It  is  anti¬ 
cipated  that  these  techniques  and  results  will  aid  in  the 
development  of  improved  models  of  the  lower  thermosphere. 
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